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Formal structure of phase-space path integrals based on different types of operator orderings is 
analysed. 



I. INTRODUCTORY REMARKS 

Section [ill of these notes was written as part of my 
discussion of phase-space path-integral methods with 
A. Polkovnikov. I wanted to translate his approach f]\ 
for myself into a familiar language. Among other things 
it is shown that if we construct a phase-space path inte- 
gral using the symmetric ordering of free-field operators, 
the path integral expresses naturally quantum averages 
of time-symmetrically ordered products of Heisenberg op- 
erators. These resuls were included, with minor changes, 
into our paper 0]. Later I realised that the resuls may 
be generalised to an arbitrary ordering of free-field op- 
erators. As is shown in section lllli for any ordering of 
free-field operators one finds the corresponding ordering 
of Heisenberg operators. In particular, using the normal 
ordering for the path integral results in Glauber-Kelly- 
Kleiner's time-normal ordering of Heisenberg operators. 
The fact that commutators of Heisenberg operators for 
different times are related to response properties of the 
system also holds irrespective of the choice of the under- 
lying operators ordering. 

At the time of writing these notes, I found it easier to 
make as much use as possible of the Weyl-based approach 
also in the case of a non-Weyl ordering. A much better 
idea would be treating all orderings on an equal basis, but 
this belongs in the future. In the meantime I decided to 
post the notes as they are. Perhaps I am not the only 
one who finds the results interesting. 



II. MULTITIME WIGNER REPRESENTATION 



so that 

A{a) ^ [A]{a) = [ ^e^°''-'^'"TTAD^P). (4) 



The alternative notation [A\ (a) for the Wigner represen- 
tation of A comes especially handj when A is an operator 
expression. Using that — Z)(— /3) it is easy to prove 

the correspondence rule for the trace, 



Tr AB = / A{a)B{a). 



(5) 



The symmetric representation of the rho-matrix p is 
called the Wigner function and often denoted W{a). 
However we follow our general pattern of using the same 
letter for the operator and its symmetric representation 
and denote the Wigner function as p{a). 

All calculations in this paper are in essence based on 
one key formula expressing the Wigner representation of 
an operator product AB, denoted as [y4i?](Q!), by the 
Wigner representations of the factors. Such relation fol- 
lows by firstly expressing the operators by their symmet- 
ric representations using eq. ([T]), and then employing eq. 
to express It is easy to verify that 



(6) 



The rest of the calculation is trivial algebra. As a result 
we obtain 



A. Preliminaries 

We firstly refresh our memory on the Wigner repre- 
sentation. We follow the paper by Agarwal and Wolf Q 
(perhaps more in spirit than to the letter). A complex 
function A{a) is the symmetric, or Wigner, representa- 
tion of the operator A if 



A 



fa 

TT 



^^^e"^*-"*'5^(a)^(/3), (1) 



where D(P) is the complex displacement operator, 
D{P)^e^^'-'^'^, Tib\l3)b(p')^Ti5^^\p- (3'). (2) 

Equation ([1]) may be obtained by using completeness of 
the set of displacement operators, 

I ^b{p)TTAD^ip), (3) 

TT 



X A{ao)B{an + a/2). (7) 



By the change of variable ao ao + a/2 we can write it 
in the alternative form 



X A(ao - c7/2)B(ao). (8) 

The usual operator-to-phase-space correspondences may 
be obtained from eqs. ^ and ([8|). For instance, from eq. 




2 



[/5a] (a) = 



p{ao){aQ + cr/2) 



Integrating by parts we move d/dag on p(cko)- In partic- 
ular this "frees" the integration over a, 



(fa 



so that the remaining integral over ao is taken trivially. 
We have thus recovered the first of the four standard cor- 
repondences characteristic of the Wigner representation, 



pa 



ap 



pa' 



a) p 



l_d_ 

l_d 
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ld_ 
2d^ 



r+2a^J 



(11) 



The other three follow with equal ease. With p[a) — > 
Aip) they hold not only for the rho-matrix but for any 
operator A. 



B. Phase-space transition amplitude 

We assume that the coincidence point for the 
Schrodinger and Heisenberg pictures is at t = to- 
The Heisenberg rho-matrix p then equals with the 
Schrodinger rho-matrix p{£) at t = Iq, 



P = /5(io)- 



(12) 



It is convenient to preserve in the notation: the Wigner 
function of the Heisenberg rho-matrix will be denoted 
p{a, to). Consider the quantum average of a Heisenberg 
operator B{t), t > to 



{B{t)) = TrB{t)p{to) 



(13) 



The expression for B{t) in terms of the evolution opera- 
tor, 



IS 



B{t) =U\t,to)B{tWM). 



(14) 



(15) 



where B{t) is B{t) in the Schrodinger picture (time- 
dependent, in general). By making use of the equations 
(P)-® we write 

{B{t))^ J ^-^2f^Bia,t)U{a,t,ao,to)pia,to), (16) 

where B{a,t) = [B{t)]{a). The phase -space transition 
amplitude U{a,t,ao,tQ) is defined by 



U{a,t,ao,to)^ J 'iJ^^e^P'-<^'P+-oK-<-Wo 

X Ttb\l3o)U\t,to)b\l3)U{t,to) (17) 
For U (a, i, ao, to) we have natural conditions 

[/(a,to,ao,to) = U{a,t,ao,to)\.f^^o ^ T^5^'^\ct - ao). 

(18) 

Unitarity of the quantum evolution is expressed in terms 
of U{a, t, ao, to) as 

(19) 



f da 

/ —^U{a,t,ao,to) = 1. 



This property is easily proven using eq. (|17p and unitar- 
ity of U{t,to). Furthermore, the group property of the 
evolution operator, [t > ti > t2) 



U{t,to)^U{t,ti)Uiti,to), 



(20) 



results in the group property of the phase-space ampli- 
tude. 



U{a,t,ao,to) 



d^ai 



U (a, t, Qfi, ti)C/(ai, ii, ao, to)- 



(21) 



Some calculations in the below are much simplified if 
we associate the group property (PTjl with the change of 
the coincidence point to ^ ^i- Namely, the Schrodinger 
picture stays put while there are many Heisenberg pic- 
tures parametrised by tp in eqs. and ([TS]). To show 
the dependence of a Heisenberg operator on the coinci- 
dence point ti we shall denote it as Bt^{t) (say). The 
absence of index signals the standard coincidence point, 

m = K{t)- 

Assume now the group property is inserted in eq. (fT^ 
We immediately recognise 



d'^ao 



C/(ai,ti,ao,to)p(ao,to) 



p(ai,ti) = [p(ti)](ai) (22) 
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as the Wigner representation of the Schrodinger rho- 
matrix p(ai,ti), which is at the same time the Wigner 
representation of the Heisenberg rho-matrix with coinci- 
dence point ti. We have then to conclude that 

J ^B{a,t)U{a,t,auh) = [Bt,{t)\{a^), (23) 

as a consequence of (f22|) and of the relation 

(i3(t))=TrSt,(t)p(ti) (24) 

expressing preservation of the averages. Using (jlSp to 
express Bt^ [t) by B{t) we find a useful relation 

B{a, t)U(a, t, ai, ti) 

TT 

= [l({h,to)B{t)U\ti,to)]{ai). (25) 

C. Phase-space path integral and the truncated 
Wigner representation 

By itself, the amplitude U{a,t, ao, to) is not associ- 
ated with a path-integral approach. In this respect it 
is similar to the coordinate-space transition amplitude, 



of which Feynman's path-integral representation is only 
one of the multitude of ways of approaching it. However, 
unlike the coordinate-space transition amplitude, a path- 
integral representation of the phase-space amplitude is a 
natural — and practical — way of looking at the phase- 
space evolution. 

While for simplicity we assumed in eq. (fT4|) that the 
Hamiltonian is time- independent, the definition (|17p does 
not make use of this assumption, nor does our notation. 
If the Schrodinger Hamiltonian H{t) is time dependent, 
in place of (fT4)) one has to use the Schrodinger equation 

im{t, to) = H{t)U{t, to), U{toM) = i- (26) 

For infinitesimally small time interval t — to = ^t the 
evolution operator is given by the approximate formula 

U{to + At, to)~i- iH{to)At/h. (27) 

One can then construct the evolution operator for finite 
time intervals as an infinite product of the infinitesimal 
evolution operators. By virtue of the group property (j2ip 
this automatically expresses the phase-space transition 
amplitude as a path integral: 



N 

U{a,t,ao,to) = lim U{a,t,aN,tN)Y\ ■ 

fe=l 



(fakU{ak,tk, ak-i,tk-i) 



(28) 



where we have assumed that the interval t — to is sliced into + 1 infinitesimal intervals At — {t — to)/{N + 1), and 
^fc = ^0 + kAt, k = 0, - ■ ■ , N are their left ends. 

Unlike the operators, all factors in ((28)) are commuting c-numbers so that the properties of the path integral can be 
fully understood from the properties of the transition amplitude over an infinitesimal time interval U {a, t + At, ao, t). 
To find it, consider the evolution of the Schrodinger rho-matrix from t to t + At: 

p{t + At) = p{t) - i[H{t)p{t) - p{t)H{t)]At/h. (29) 

By making use of eqs. ([7]) and ([8|) we can turn this into a relation between the corresponding Wigner functions: 



p{a,t + At) = 



(Q-an)(T*-(Q-cto)*o- 



lAt r 

p{ao,t){l - — H{ao - <7/2,t) - H{ao + <j/2,t) \, (30) 



where H{a, t) is the Wigner representation of the Hamitonian. The zeroth-order term was dragged under the integral 
using eq. (fTU|) . To recover U{a, t + At, ao, t) from this formula we note that it follows from that 



—^U{a, t + At, ao, t)p{ao, t). 



This way. 



iAt r 



U {a, t + At, ao,t) — / exp < (a — ao)a* — (a — ao)*<7 — H{ao — a/2,t) — H{ao + cr/2, t) 



(31) 



(32) 



In this formula, we have also replaced the expression in the curly brackets in (|30p by the exponent. This is justified 
for an infinitesimally small At. 
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Following [l| we separate in (|32p the classical evolution from quantum fluctuations. Formally, we write 

H{ao + a/2, t) - H{ao - a/2, t) = a/*(ao, t) + a*f{ao, t) + h<^^\ao, a, t), (33) 

where 

dH{ao,t) dH{ao,t) 

This way, 

+ Ai,Q;o,i) = J ^ exp |(q! - ao + ifAt/h)a* - {a - ao +ifAt/hya + ^/i(^)(ao, cr, t)| . (35) 

For polynomial Hamitonians /i*-^-* (ag, cr, i) is a polynomial where all terms are at least cubic in a, a*. Nonzero 
h'-^^ao, a, t) can only occur for nonlinear interactions. If for some physical reason can be neglected, then 

t/(a,< + Ai,ao,i) = Ti5'^'^\a - + if At/h). (36) 

In this case the motion in phase space is along the trajectories obeying the equation 

iha = f{a,t). (37) 

This motion is deterministic (nonstochastic) and in this sense classical Nonzero ft.'^^-' introduces quantum fluctua- 
tions [8] which may be accounted for exactly _4] or perturbatively 



D. Two-time averages 



Consider the quantum average (^o < ^1,^2) 



(38) 



where ^{t) and y{t) are Heisenberg operators. The evolution operator relates them to their Schrodinger representa- 
tions X, Y, cf. eq. Rather than distinguishing the cases ti > t2 and ti < t2, we assume that ti < t2 and consider 
two distinct averages, {X{ti)y{t2)) and {y{t2)X{ti)'). Simple formulae may be found for the practically important 
case X{t) = A{t),A''{t), where 

A{t) ^U''{t,to)aU{t,to), A\t)^U\tM)a)U{t,tQ), (39) 

are the Heisenberg field operators. Consider, for example, the average {y{t2)A{ti)) . Moving the coincidence point to 
t ~ ti and using eq. we have 



{yit2)A{h)) =TTytAt2)ap{h) = 



(Pai(Pa2 



Y{a2)U{a2,t2,ai,ti) [ap{ti)] (ai). 



The phase-space correspondences ([TT|) allow us to write 



[apih)] (ai) = (ai -|- no- 
where p{ai,ti) is given by eq. (|22p . As a result of these manipulations we obtain 



{y{t2)A{ti)) 



(Pa2(Pai(Pao 



Yia2) 



(ai - ^-^)u{a2,t2,ai,ti) 



U{ai,ti,ao,to)p{ao,to) 



(40) 



(41) 



(42) 



Integration by parts was used to move the derivative to U{a2, ^2, cti, ii); square brackets emphasize that the differen- 
tiation does not apply to [/(ai, <i, ao, io)- Similar considerations yield 



{A{t^)y{t2)) - 

{A\h)y{t2)) 

{y{t2)A\h)) 



(Pa2(Paid?ao 
(Pa2(Pai(PaQ 



Y{a2) 



Y{a2) 



Y{a2) 



1 d 
29^^[ 

1 d 



f/(a2,i2,ai,ii) 



/ * . 1 5 



2dai 



U{a2,t2,ai,ti) 



U(ai,ti,ao,to)p{ao,to), (43) 
U{ai,ti,ao,to)p{ao,to), (44) 
U{ai,ti,ao,to)p{ao,to). (45) 



We remind the reader that eqs. hold if io < < ^2- Note that the latest operator in them remains arbitrary. 
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E. Commuting Heisenberg operators as the quantum response problem 

Relations ([43l) - (|45)) are exact and not associated with the path-integral representation of the phase-space amplitude. 
However their most natural interpretation is in terms of path-integral averages. Adding (|43p and (I42p together we 
find 

^{yit2)A{h) + A{ti)y{t2)) = J ^-^^^-^^^-^Y{a2)aiU{a2,t2,ai,ti)U{ai,ti,ao,to)piao,to). (46) 

This is nothing but the path-integral average of Y {a{t2))a{ti) 9], where a{t) is an individual path. Hence 

l:{yit2)A{h) + A{h)yit2)) = Y{a{t2))a{h), 

T ^ (47) 

^{y{t2)AHti)+AHh)y{t2)) = Y{a{t2))a*{h). 

where the bar denotes the (quasi-statistical) averaging over the paths. The second line here follows by adding together 
(|44)) and p5|) . These relations suggest that the path-integral averages in the Wigner representation correspond to 
quantum averages of symmetrised products of Heisenberg field operators. Natural as it looks, this result in fact holds 
only for two-time averages. For three and more times, a new type of operator ordering is found, see section fllF I 

It is easy to see that derivatives in eqs. (|^ - (|^ correspond to commutators of the Heisenberg operators. For 
example, subtracting ([44]) from ([45]) yields 

/f<ju \ aU^ M\ f d?OL2(Paid?ao , .dU{a2,t2, ai,ti) 

{[y{t2),A^{ti)\) ^ / -3 Y{a2) U{auti,ao,to)p{aQ,to), (48) 



As was suggested by Polkovnikov [1], quantities like ((48|) 
may be calculated numerically within a path-integral ap- 
proach by introducing "quantum jumps" in the trajecto- 
ries. In practice, this means splitting a trajectory a,t t — 
ti by introducing small complex shifts a{ti) — > a{ti)~\-5a, 
in order to calculate the derivative numerically. For a 
complex derivative, two independent shifts are needed, 
hence for t > ti three trajectories must be run in par- 
allel. The numerical cost of such "quantum jumps" is 
obviously not prohibitive. 

A similar method was employed in Q to calculate sym- 
metrised averages of Heisenberg operators using the so- 
called positive-P representation. Path-integral averages 
within the latter correspond to time- normally ordered op- 
erator averages. The commutators in [5] were expressed 
by, firstly, using Kubo's famous formula for the linear 
response function [Q] to associate the commutator with 
response of the quantum system, and, secondly, connect- 
ing the latter to a response problem in phase space. It is 
easy to see that such response interpretation also holds 
for the derivatives in eqs. ([42l) - (|45)) . Indeed, assume that 
the Hamiltonian of the system has been complemented 
by external c-number sources. 



H{t) H{t) - s{t)a) - s*{t)a. 



(49) 



In phase-space terms, this only modifies the regular evo- 
lution, 



/(a,t) 



r{a,t)-s 



sit), 



it), 



(50) 



cf. eqs. Shifting a trajectory by Sa at t ~ ti 

requires an instantaneous source 



s{t) = —ih5ad{t — ti). 



(51) 



We thus have a simple correspondence between sources 
and "quantum jumps:" (cf. also endnote @) 



I d 
h dai 



5s* [ti 



i__d_ 

hda\ 



The RHS of (gS]) is then recognised to be, up to a factor, 
the linear response function of the quantum system. We 
can therefore rewrite (1^5]) as 



5s{ti) 



= j^{[y{t2),AHh)\ 



t2 > U 



(53) 



This is nothing but Kubo's formula for the linear re- 
sponse function. Hence equation (j48p can be derived 
starting from Kubo's formula and then using the cor- 
respondence between the sources and "quantum jumps." 
A general discussion of the quantum response problem 
in the Wigner representation requires advanced formal 
tools and will be presented elsewhere. 



F. Multi-time averages and time-symmetric 
operator ordering 

Consider the path-integral average {ti <t2 < ■ ■ -1^) 



a{ti)a{t2) ■ ■ ■ a{tN) = {TwA(ti)A{t2) ■ ■ ■ A{tN)), (54) 
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where the Heisenberg field operators are given by ([39)) . 
The LHS of eq. is defined in phase-space terms as 



a{ti)a{t2) ■ ■ ■ a{tN) ^ f ""'^ '^li+i ^ "^ ^("0,^0) ctkU{ak,tk,ak-i,tk-i)- (55) 

fc=i 

On the RHS of (|54p we redefine it in the Hilbert-space terms as an average of a time- symmetric product of the field 
operators. 

To obtain an independent characterisation of the time-symmetric ordering we use the same trick that lead to eqs. 
(|42)) - ([45|) . but apply it "in reverse." We start from defining the quantity 

/(Pa (Pa ^ 
^ — —a2Y[cfkU{ak,tk,ak-i,tk^i). (56) 

fe=3 



The times ts, - ■ ■ ,tpf here should be considered as parameters. Following the pattern of eqs. ([23)1 and (|25|) we then 
write 

(P 

-^B{a2,t2)U{a2,t2,ai,ti) = [S*, (^2)] (ai) = [W(ii, to)^(i2)^^(ii, ^o)] (^i)- (57) 
Furthermore, by making use of eq. (|22p and of the phase-space correspondences, we have 

/d^ (y. 1 
—^U{ai,ti,ao,tf))p{aQ,tQ) = aip{ai,ti) = - [ap{ti) + p{ti)a\{ai) (58) 

The remaining integration over ai in (|55p is then performed straghtaway and we find 

{TwA{ti)A{t2)---A{tN)) - iTrW(ti,to)S(t2)Wt(ti,io)[a/5(ti) + p(ti)a] = ^({^(ta), i(ti)}), (59) 
where the curly brackets stand for the anticommutator, 

{X,Y]=XY + YX. (60) 

The last equality in ((59l) follows from 

p{ti) = U{ti,to)p{to)UHti,to)- (61) 

Finally, the operator 6(^2) is found if we take eq. ([57]) at ti — to and compare it with eq. (j55|) . It is then easy to see 
that 

B{t2) = TwA{t2)---A{tN). (62) 

Equation ([55]) is thus a recursive relation for the time-ordered operator products. 

If we replace in ([Ml) a{ti) a*{ti) and A{ti) A''(ti), equation ([59]) will also hold with A{ti) A^{ti). 
Furthermore, the actual nature of the factors a(t2) ■ ■ ■ a(t]\r) is irrevant. Any subset of them may be complex- 
conjugated, provided the corresponding operators under the Tiy-ordering are Hermitian-conjugated. As a result, we 
arrive at the following recursive definition of the time-symmetric ordering: 

Twi = i, 
TwA{t) = A{t), TwAHt) = A^t), 

TwA{t)V>t^l{A{t),r>t}, ^^^^ 
TwAHt)Vyt^l{A\t),Pyt}, 
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where P>t is a time-symmetric product of operators with all time arguments larger than t. 

The number of different terms in a time-symmetric product of iV factors is 2^~^. For N > 2, this is clearly different 
from iV! permutations comprising a symmetrised product of the same factors. The time-symmetric and symmetrised 
products are thus distinct for all N > 2. For example, for = 3 and ti < t2 < we have 

TwA{t,)A(t2)A{t3) = i [A{h)A{h)A{h) + A{t2)A{t3)A{h) + A{ti)A{h)A{t2) + A{t3)A{t2)A{h)] . (64) 

The symmetrised product of the same factors should also include A{t2)A{ti)A{t3) and A{t3)A{ti)A{t2)- 

The most important properties of the time-symmetric products are the two: these products are continuous at 
coinciding time arguments, and for free-field operators they turn into the conventional symmetric (Weyl) ordered 
products. A proof of these properties, as well as a more detailed discussion of the time symmetric ordering, including 
its relation to Schwinger's closed-time-loop formalism, will be presented elsewhere. 

III. MULTITIME REPRESENTATIONS BASED ON OTHER ORDERINGS 

A. Definitions 

Non-symmetrically ordered representations are introduced by replacing eq. (jj]) by 

A.(a) = [i]^(a)= I ^?!^e'3"*-/''"+^l/'lV2Tri^t(^)^exp(-^^^)A(a). (65) 

Equation ([5]) then becomes 

J - f d^c( 
TtAB = — 

r 



exp 



9^ 



2 dada* 



A{a) 



exp 



2 dada* 



B{a) 



-As{a)B.s{a). 



(66) 



The normal, Weyl and anti-normal representations of the 
p-matrix, corresponding to s = —1,0,1, are commonly 
known as Q, W and P functions, respectively. They 
are also called quasi-distributions, because, for example, 
the P-function gives a classically-looking expression for 
averages of normally-ordered operator expressions, 

/i2 
-^P{a)F{a,a*), (67) 

and similarly for the Q-function with the antinormal op- 
erator ordering. Thus normal representation of field op- 
erators is associated with antinormally ordered p-matrix, 
and vice versa. To avoid lengthy repeated specifications, 
we introduce the term s -representation for a phase-space 
picture where the field operators are s-ordered and the 



/9-matrix is (— s)-ordered. This agrees fully with the use 
of this term in quantum optics: normal representation 
employs P-function while antinormal representation em- 
ploys Q-function, and so on. 

B. Pliase-space transition amplitude 

Manipulating similar to yields in the s- 

representation 

/ri/ \\ f <PaQ(Pa 

\o{t)) = / ^ Bs(a,t)Us[a,t,ao,to)p-s{a,to), 

(68) 

where Bs{a,t) = [P(t)]^(Q!) and 



Us{a,t,ao,to) = exp 



92 



s / d 

2 \dada* ' daedal 



U{a, t, ao, io) 



is the phase-space transition amplitude in the s-representation. It is easy to see that the properties (IT51) - (PT|) also 
hold for Us, while eqs. ([22])- ([251) become 



d'^ao 



Us{ai,ti,ao,to)p^siao,to) = P-s(ai,ti) = [p{ti)] _^{ai) 



(70) 
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and 

-^B,{a,t)Usia,t,a,,h) = [^t, (i)] ,(ai) = [l(ih,to)B{t)U^ ih,to)] ^{a,). (71) 



C. Path-integral representation of the amplitude 

While the general structural properties of the path integral in phase space do not change with the ordering, the 
path integral itself is certainly ordering-specific. Consider the infinitesimal amplitude U{a,t + At,ao,t). Combining 
eqs. ([32]) and ((69l) we have 



Us{a, t + At, aQ,t) — exp 



92 



92 



-f- 

2 \dada* daoda^ 



■ exp [(a — oiQ)a* — (a — Q;o)*cr 



X 1- 



Hiao-(7/2,t)-H{aa + c7/2,t) }. (72) 



A straightforward calculation in Appendix |A] yields 



Us{a,t + At,ao,t) 



exp s (a — ao)""* — {a — q:q)*(7 



iAt 



Hsiao- 



(1 + s)a , (1 - s)a* 



, i -Hsiao 



(1 - s)a 



"0 



(l + s)a* 



(73) 



We had to amend our notation, 



Hs{a,t) ^ Hs{a,a\t), 



(74) 



to accommodate for asymmetry between a and a* introduced by a non-symmetric ordering. Separating the deter- 
ministic part of the motion in phase space we get in place of eq. (j35p 



Us{a, t + At, ao, t) 



exp i (a — ao + ifsAt/h)a* — (a — ao + ifsAt/h)*a + ——hs{ao, a, t) 



(75) 



where 



dHs{ao,t) dHs{ao,t) 



(76) 



and 



(1 - .s)cT , (1 + ^)^* , 
"0 H , "o ^ ' * 



- Hs\aO , Qfg 



,t) -(7/;(«o,t)-^*/.(«o,t). (77) 



Expanding this into power series we get for the practically important cases of the normal, antinormal, and Wigner 
representations 



/ii(ao, cr,t) = ^ 



fe>2 



/i_i(ao,cr,t) = ^ 



fe>2 



fc! ' 
fc! ' 



/io(ao,a',t) = ^ 



/c-t-m>3,odd 



2'=+™-ifc!TO! ' 



(^) 

(Q) 



(78) 
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where 



[/^™)(ao,t)]*=/i' 



(79) 



,k) 



{aQ,t). 



Unlike in the symmetric case, for an arbitrary s the noise contribution starts from the second-order noise. Interest- 
ingly, the deterministic path is also ordering-specific. The classical motion is only recovered in the limit h ^ 0, when 
both the noise contribution, and the reordering correction in which Hs{ao,t) differ for different s, are negligible. 



D. Quantum averages 



In the s-representation the phase-space correspondences read 

.s + 1 d 

pa 
cip 
pa 



a) p 



(° 

'-(' 
(' 



2 


da*) 


s - 1 


' ) 


2 


da*) 


s — 


1 d \ 


2 


da) 


s + 


1 d \ 



2 d 



)p-s{a). 



(80) 



These formulae may be derived from the definition (|65p . In fact, the coefficients at the derivatives may be deduced from 
eqs. if noticing that, for s = 1 (s = —1), the p-matrix is antinormally (normally) ordered, hence correspondences 
for dp and pa) {pa, and a^ p) should be without derivatives. Furthermore, the way eqs. ((80|) are written emphasizes 
that the "nonsymmetric correction" in these always looks as — s/2 times a derivative. The route from the standard 
correspondences to eqs. (j42 |) -(|45 |) involves one integration by parts, hence all "nonsymmetric corrections" to (j42|) -((45 |) 
must look as s/2 times a derivative. This allows us immediately to write 



{yit2)A{t,)) 
{Mtim2)) 

{y{t2)A\h)) 



cPa2cPaid?aQ 



tPa2(Pai(fao 



Ysia2) 



Ys{a2) 



-1 d 



ai 



ai 



(Pa2(Pai(PaQ 



(Pa2(Pai(PaQ 



Ys{a2) 



Ysia2) 



2 


dal) 


s + 


1 d X 


2 


dal J 


s — 


1 d ^ 


2 


dai ) 


s + 


1 d \ 



Us{a2,t2,ai,ti) 



Usia2,t2,ai,ti) 



Usia2,t2,ai,ti) 



Usia2,t2,ai,ti) 



Us{ai,ti,ao,to)p-s{ao,to). (81) 



C/s(ai,ii,ao,to)p-s(ao,io), (82) 



t/s(Q;i,ti,Q;o,io)P-s(ao,io), (83) 



Usiai,ti,ao,to)p_sio:o,to). (84) 



' 2 dai. 

E. "Quantum jumps", response, and the time-s-ordering of Heisenberg operators 

The property that adding external sources to the Hamiltonian is equivalent to adding the same sources to the 
equation for the deterministic path in phase-space does not depend on the operator ordering, hence eqs. (|48p - (|53p are 
valid for any s. This "invariance of the response viewpoint" holds the fact notwithstanding, that the path-integral 
averages themselves correspond to different types of quantum averages for different values of s. On top of purely 
structural properties of the path integral, the reasoning in section Hi Fl is based on the relations 



i [ap + pd] (a) = ap{a), ^ [d^ p + pd'^] {a) = a*p{a), 



(85) 



10 



characteristic of the Weyl ordering. The corresponding relations in the s-representation read 



1 



1 - s 



-ap - 



-pa 



(a) = 



1 



1 



-a' p 



-pa' 



{a) = a* p-s{a). 



(86) 



We now define the time- s- ordered product of the Heisenberg field operators by generaling eq. (|54p to the case of 
s-representation, 



y{h)a{t2) ■ ■ ■ ait^r' = {%A{h)A{t2) ■ • •i(tjv)). (87) 
The path-integral average on the LHS here is given by eq. ([55)) with U ^ Ug and p — > p-g', full meaning to it is then 



assigned by eq. ([73]) . The first relation in IIIFI that has to be modified beyond these trivial replacements is eq. ([58 



the last equality in it must follow ([86)1 instead of 



As a result the recursion relation ([55)1 changes to 



1 + 



B^'Ht2)Aih) + ^Aih)B'^'Ht2 



with 



B^'Ht2)^%A{t2)---A{tN). 



(89) 



Similar reasoning applies with A{ti) A''{ti), except 
that with this replacement one should also swap the 
weight factors in ([55)1 . We thus arrive at the recursive 
definition of the time-s-ordered product: 



T,i = a, 
%Ait) = Ait), TsAHt) = AHt), 



%Ait)rli 
%A\t)r^^l : 



1 



2 ' ' 2 
1 + s 



(90) 



" is) 

where V^l is a time-s-ordered product of operators with 
all time arguments larger than t. 
Properties: 

• For s = 1 and s = 0: time-normal and time- 
symmetric products. 

• For other s — new orderings; e.g., time-antinormal 
for s = — 1. 

• Hermitian-conjugate of a 7^-product is a 
product. 



±lA\t)V^^l + ^-^V^lA\t), 



Continuity? (hopefully) 



J 



APPENDIX A: INFINITESIMAL TRANSITION AMPLITUDE IN THE S-REPRESENTATION 



Here we outline the reasoning leading from eq. ()72p to ([75)1 . We assume that the differential operator in ([7^ may 
be drawn under the integral and applied directly to the two factors comprising the integrand. Applying exp (| g^g^* ) 
is straightforward and results in a factor exp ( — ^^^^ ). Applying exp ( — § daoda" ) ^^^"^ easy if using the formula 
(with Fi , F2 and being arbitrary complex functions) 



d d 



d 



d d 



d 



^1 7^' TTT K2(ao)F3(ao) =Fi — + —, — + F2(ao)F3(a[,)|^,_^, 
which expresses the rule of product differentiation. To start with, we write 



exp 



s / d (J \r a 



d 



d 



2 V dao da^ / V da^ da^ 
= exp 

Then, firstly, 



92 



2 \daoda'Q da^daQ 



exp 



92 



2 daoda} 



exp 



2 da'f^da'^ 



exp 



92 



2 danda, 



— J exp [(a — ao)(7* — (a — ao)*"'] = exp [(a — aa)<7* — (a — ao)*^ + ■s|o'|^/2] . 
0/ 



(Al) 



(A2) 



(A3) 
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This cancels exp ( — ^^y-) emerging from exp (| q^q^, ) ■ Secondly, from eq. (I65p 



exp 



2 da'^^da'* 



1 - 



- cr/2, <) - + a/2, t) 



«At r 



iJ.K - a/2,t) - i/,(a^ + a/2,t) . (A4) 



Thirdly, by using the commutational relations for derivatives, 
d 



we have 



_d 

5a, 



exp 



exp [(a — ao)(T* — (a — ao)*cr] — exp [(a — ao)(7* — (a — ao)*o'] 



9an 



- exp [{a - ao)(T* ~ (a - aQ)*a] = exp [{a - ao)(J* - {a ~ ao)*a] ( ^— ^ + a 



X 1 



exp [(a — ao)o'* — (a — ao)*o'] exp 



exp [(a — ao)cr* — (a — ao)*cr] 
iAt r 

ff.K - (7/2, t) - H,{a'o + a/2,t) 
sa* d sa d 



2 da^ 2 



exp [(a — ao)""* — (a — ao)*o'] <{ 1 — 



/ (l + g)g , (l-g)fT* , 
rls[oio , an ^ ^ 



-i/,(ao + ^ 7r^,a* + ^ :r^,t 



In this formula, we had to amend our notation. 



Hsia,t) ^ Hs{a,a*,t), 



(A5) 



(A6) 



(A7) 



to accommodate for asymmetry between a and a* introduced by a non-symmetric ordering. We have thus indeed 
recovered the integrand of eq. ([73| . 
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